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We present deep polarimetric observations at 1420 MHz of the European 
Large Area ISO Survey North 1 region (ELAIS Nl) as part of the Dominion 
Radio Astrophysical Observatory Planck Deep Fields project. By combining 
closely spaced apertures synthesis fields, we image a region of 7.43 square degrees 
to a maximum sensitivity in Stokes Q and U of 78 /iJy beam~^, and detect 786 
compact sources in Stokes /. Of these, 83 exhibit polarized emission. We find 
that the differential source counts (log A^ — log p) for polarized sources are nearly 
constant down to p > 500 yuJy, and that these faint polarized radio sources are 
more highly polarized than the strong source population. The median fractional 
polarization is 4.8±0.7% for polarized sources with Stokes / flux density between 
1 and 30 mJy; approximately three times larger than sources with / > 100 
mJy. The majority of the polarized sources have been identified with galaxies in 
the Spitzer Wide Area Infrared Extragalactic Survey (SWIRE) image of ELAIS 
Nl. Most of the galaxies occupy regions in the IRAC 5.8/3.6 fim vs. 8.0/4.5 /xm 
color-color diagram associated with dusty AGNs, or with ellipticals with an aging 
stellar population. A few host galaxies have colors that suggests significant PAH 
emission in the near-infrared. A small fraction, 12%, of the polarized sources are 
not detected in the SWIRE data. None of the polarized sources in our sample 
appears to be associated with an actively star-forming galaxy. 

Subject headings: polarization — techniques: polarimetric — radio continuum : 
galaxies — galaxies : evolution — individual (ELAIS Nl) 



1. Introduction 



Observation of polarized radiation at radio wavelengths is one of the prime means to 
study the roles of magnetic fields in astrophysics, through synchrotron emissivity which 
samples magnetic fields in relativistic plasmas, and through Faraday Rotation, created by 
radiation propagation through magnetized thermal plasmas. Understanding the origin and 
evolution of magnetic fields is a key science goal of the The Square Kilometre Array, a 
next- generation radio telescope under development by the international community (ISchilizzi 
20041 ). The primary observational tool to study the magnetic universe will be an SKA all-sky 
Rotation Measure survey of background ra dio sources and diffuse Galactic emission down to 
polarized flux density levels of ~ 0.1 yuJy (IBeck &: Gaenslerll2004l ). However, w hile modern 
source counts approach fl ux density sensitivities of ~ 10 /xJy in total intensity (IWindhorst 



20031 : iHopkins et al.ll2003l ). very little is known about the polarization properties of the faint 
radio source population. The most extensive analyses of polarization of compact extragalac- 
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tic so urces were carried ou t by lMesa et al.l (120021 ) and lTucci et al.l (120041 ). who used the NVSS 
data (ICondon et al.l Il998l ) to derive statistical polari z ation properties for ~30,000 sources 
with S'i,4GHz > 100 mJy. Similarly, iBeck &: Gaenslerl (120041 ) used NVSS sources with total 
flux density greater than 80 mJy to extrapolate polarized source counts to yuJy levels. 



Mesa et al.l (120021 ) found that the mean fractional polarization of radio sources in the 



NVSS brighter than 80 mJy was anti-correlated with flux density , especially for steep 



spectrum radio sources (a < —0.5, for Sy ~ z/"). iTucci et al.l (120041 ) conflrmed this result 
for the median of the fractional polarization for steep-spectrum sources only (87% of their 
sample), but foun d no signiflcant trend for flat spectrum sources (13% of their sample). 
Tucci et al.l ( 120041 ) also noted that the flat shape of the polarized source counts indicates a 
dependence of the fractional polarization on flux density. 

Radio sources with 5'i.4ghz > 100 mJy are predominantly associated with flat or steep 
spectrum Active Galactic Nuclei. Star forming galaxies begin to be a signiflcant fraction of 
the population at flux densities less than a few mJy (see e.g. Hopkins 2000, Windhorst 2003). 
However, there is still ongoing debate abou t the fraction of the sub-mJy radio sources which 
is radio-quiet active galactic nuclei (AGN) (iGruppioni et al.lll999l : ISimpson et al.ll2006l ). It is 
thus uncertain that polarization properties derived from the strong radio source population 
are applicable to the sub-mJy radio sources. 

We have begun a project called the DRAO Planck Deep Fields to explore the high 
latitude sky at high sensitivity in polarized radio continuum and in atomic hydrogen emission. 
The project uses the DRAO Synthesis Telescope at 1.4 GHz to create deep images of two 
flelds, one with a very low column of foreground material, the ELAIS Nl region (/, h) = 
(84°, +45°), and a larger region of highly structured infrared cirrus emission at {l,b) = 
(135°, +40°). This paper reports initial results from the flrst 30% of o bservations of ELA IS 
Nl. The ELAIS Nl (European Large Area ISO Survey North 1) fleld Joiiver et al.lboooh . is 
an area of approximately 2 square degree chosen for a mid-infrared survey of distant galaxies 
with ISO. The fleld was selected to minimize confusion with Galactic cirrus and the zodiacal 
background, being one of the lines of sight to the extragalactic sky with minimum IRAS 100 
/im emission. A larger area that includes ELA IS Nl was later observe d by the Spitzer Wide 
Area Infrared Extragalactic Survey (SWIRE) (iLondsdale et al.ll2003l ). The large amount of 
archival data from these extragalactic surveys make this region ideal for studies of the faint 
polarized radio source population. 
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2. Observations and Data Processing 
2.1. Synthesis Observations 



The DRAG Synthesis Telescope (DRAG ST) is described in detail in iLandecker et al. 



(|2000| ). The telescope is a seven-element east-west array of 9-m diameter antennas. Three 
antennas are moved to provide complete sampling of the UV plane from the shortest baseline 
(12.86 m) to the longest baseline (617.18 m) after a full synthesis of 12 times 12 hours. The 
array has a primary beam size 107'2 (FWHM) at 1420 MHz, which makes it an effective 
instrument for wide-field surveys. The first sidelobe of the synthesized beam is at the 3% 
level, and side-lobes farther from the main lobe of the beam are less than 0.5%. The first 
grating ring of the synthesized beam appears at 2?8 from the main lobe at 1420 MHz, which 
is outside the field of view defined by the 10% sensitivity level of the primary beam. 

The telescope observes simultaneously the H I 21 cm line and continuum at 408 MHz 
and full polarimetry in four 7.5 MHz wide frequency bands centered around 1420 MHz. The 
antennas have prime focus feeds and at 1420 MHz receive both right-hand (R) and left-hand 
(L) circularly polarized radiation. The observations and data processing techniques used in 
this paper follow those that have been used to obtain high-fidelity wide-field polarimetri c 



images of the Galactic plane with the DRAG ST (ITaylor et al.ll2003l : ILandecker et al.ll2007l ). 



The system temperature of the telescope as described by ILandecker et al.l (120001 ) was 
60 K, leading to an rms noise in a synthesized image of a single field of 280 /xJy beam~^ at 
field center. However, starting in 2003, the sensitivity was enhanced by a series of improve- 
ments, completed by the time the current observations began in 2004. These improvements 
comprised installation of new low-noise amplifiers, modifications to telescope structures to 
reduce ground noise, and installation of shielding fences to further block ground radiation 
from entering the aperture. The system temperature of the seven individual antennas now 
spans the range 35 K to 60 K, and the overall system temperature is ~ 45 K, leading to a 
field-center rms noise of 210 /zJy beam~^ (53 sin 6 mK) at 1420 MHz. A higher sensitivity and 
larger field of view are obtained by creating a mosaic of overlapping fields. The most uniform 
sensitivity across a mosaic is obtained if the fields are centered on a hexagonal grid. The 
field centre separations for the survey presented here is 22', which is 20.5% of the FWHM 
diameter of the primary beam. 

Observations for the DRAG ELAIS Nl survey began in August 2004 and will continue 
to create a final mosaic of 30 fields. This paper presents 21-cm continuum polarimetry of the 
first 10 fields. Figured] shows the location of the 10-field mosaic in relation to the ISO ELAIS 
Nl field, and the SWIRE survey of ELAIS Nl. The theoretical maximum 1-a sensitivity in 
our mosaic of 10 fields is ~ 80 /iJy. 



- 5 - 



2.2. Polarization Data Processing and Calibration 

Complex gains for the center of each field were determined by observing the unresolved 
and unpolarized sources 3C 147 and 3C 295 between 12-hour observing runs. The absolute 
polarization angle was calibrated by observing the highly linearly polarized source 3C286 
once every 4 days. 

The polarization images do not contain much flux, making self-calibration ineffective. 
Therefore R and L gain solutions derived from self-calibration of the Stokes / images were 
applied to the polarization data as well. The resulting visibility data sets for each field were 
then corrected for the effects of instrumental polarization across the field of view (which leads 
to leakage of Stokes / power into Q and U). Instrumental polarization is a complex quantity, 
with phase and amplitude terms. These were measured for the seven individual antennas 
with a holographic technique, using the unpolarized source 3C 295, on a 15' grid across the 
primary beam. After interpolation, these measurements were used to predict conversion of / 
into Q and U at any point in the beam based on CLEAN component models from processing 
the / image for a field. Residual errors in instrumental polarization for an individual field 
are estimated at 0.25% for the field center, growing to 1% at a distance of 75'. With a field 
separation in the mosaic image of 22' the instrumental polarization in the central 3.6 square 
degrees of the mosaic remains less than 0.5%. Towards the edge of the mosaic instrumental 
terms may be as large as 1%. 

After these initial procedures there are usually still confusing arc-like structures left in 
the images, centered on bright sources inside and outside the primary beam. These are the 
result of residual complex gain errors at large distance from the field center. Effects from 
these are removed from the visibilities using a proc edure called modcal, which is in principle 



a direction-dependent self-calibration (IWillislll999l ). 



The antenna sidelobes are highly polarized, and sources outside the primary beam can 
produce strong spurious polarized signals. The Sun is seen in the sidelobes whenever it is 
above the horizon, but its effects are usually confined to the shorter baselines because of its 
large extent. The effects of the Sun were removed by making images centered on the Sun's 
position and removing the response from the visibilities. Terrestrial interference, which is 
always polarized, is another source of spurious polarization. Interference effects in images 
appear as a spurious source concentrated around the North Celestial Pole, and can be largely 
eliminated by making an image at the Pole and correcting the visibilities. Radiation from 
the ground appears to be polarized, and can appear as a correlated signal in visibilities 
corresponding to short baselines; this effect is more difficult to remove. In some particularly 
bad cases the data for affected interferometer spacings were simply flagged and removed. 
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2.3. The Images 

Figure [2] shows the deep 21 cm continuum images of the ELAIS Nl area in Stokes /, 
Q, and U . The images are centered on ^2000 = lO*^!!™, ^2000 = +55° and cover an area of 
7.43 square degrees (the area within the thick gray hne in Figured]). The angular resolution 
varies with the declination over the images and is given by ba x bg = 49" x 49"cosec5. At 
the mosaic center the resolution is 49" x 59". The noise near the center of the mosaic (white 
square in Figure [2]) is measured to be 78 /iJy in Q, and U (see § 13. ip . The Stokes I image 
is not limited by confusion, but the noise is slightly higher at 85 /iJy near the center of the 
mosaic, probably because of a contribution from faint sources. 

The final images are virtually free of artifacts, so the sensitivity is limited by the noise. 
The rms brightness sensitivity is 17.4 mK. The dynamic range near the center of the mosaic 
is more than 3000:1 in Stokes /. With very few exceptions, sources in Figure [2] appear as 
compact (nearly) unresolved sources, that can be characterized by their peak intensity and 
a single polarization angle. This is consistent wi th expectations fr om the angular size - flux 



density relation for extragalactic radio sources (jWindhorstI 120031 ) : the median angular size 
of a IJy source is ~ 10". 



3. Compact Polarized Sources 
3.1. Source Detection 



Flux densities and positions of all sources in the pilot deep field images were measured 
with a source extraction algorithm that fits a two-dimensional Gaussian and a background 
level to each source. The mosaic images were multiplied by the primary beam response 
function of the mosaic to obtain an image with a uniform noise level, equal to the noise level 
at the center of the image. This operation retains the correct signal-to-noise ratio for each 
source, but the resulting uniform noise level greatly facilitates automated source extraction. 
The inverse primary beam correction is applied to the measured flux densities of sources 
from the uniform noise map to transform back to true flux density. The rms noise level in 
the uniform noise images was measured by fitting a Gaussian function to the distribution of 
amplitudes in the image. The result is shown in Figure [31 The noise distribution in both 
the Q and U images is well fitted by a Gaussian with dispersion a = 78 yuJy. 

The polarized flux density image {p = \/Q^^\-TP) was searched for polarized sources. 
For Gaussian noise in the Q and U images with rms amplitude a, the statistical distribution 
of the noise in a measurement of p for a source with an intrinsic polarized flux density of Po 
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is a Rice distribution (IRi call 9451 : IVinokurlllQGSl : ISimmons &: Stewartlll985l ) 



= ^-e-'^U?E^) , (1) 

Here Jo is the modified Bessel function of the first kind. For po = 0, the distribution reduces 
to a Rayleigh distribution, 

f{p\0) = ^e-^ , (2) 
a 

which gives the probability distribution of pixel amplitudes in the p image in the absence of 
polarized emission. The noise in the p image has a non-zero mean and has higher probability 
of positive peaks above a given detection threshold than Gaussian noise. We searched the p 
map down to a level of 4.55a, which has an equivalent probability for false positive signals to 
the 4(j level for a Gaussian distribution. The measured polarized flux density p was corrected 
for noise bias to obtain an estimate of po through the relat ion vi = jp' — a"^, which is a good 



approximation if the signal to noise ratio is larger than 4 (ISimmons fc Stewartlll985l ). 



The 83 sources detected are listed in Table [H which gives the position of each source, 
the integrated total flux density, noise bias-corrected peak polarized intensity, polarization 
position angle, fractional polarization defined as the ratio of the bias-corrected peak polarized 
intensity to the peak total intensity, and the spectral index of the total flux density between 



325 M Hz and 1420 MHz if the source appears in the WENSS catalogue (iRengelink et al. 



19971 ). To ensure possible spurious polarized sources due to instrumental polarization are 
not included, we conservatively remove sources with observed fractional polarization, H, less 
than 1%. Only two sources were removed from the sample for this reason. 



3.2. Distribution of Fractional Polarization 



The intrinsic fractional polarization Ho of radio sources provides astrophysical informa- 
tion about the nature of the polarized sources. However, the observed fractional polarization 
H is sensitive to the noise in p and in /. In addition to p being a biased estimator of Po, 
the error distribution of the ratio p/I has strong non-Gaussian wings, so H is not a very 
accurate estimate of Ho even for relatively high signal-to-noise sources. This is illustrated 
in Figure HI which shows the relation between H and Stokes / flux density for an artificial 
sample of sources, all with Ho = 5% and Gaussian noise added with equal amplitudes in 
Stokes /, Q, and U . Only those sources with a p flux density more than 5a are shown. Error 
bars represent standard error propagation in H from the errors in p and I. The high values 
of H at signal-to-noise ratio less than ~100 is a result of the detection threshold in p and 
the non-Gaussian error statistics of H. The effect is much larger than the polarization noise 
bias alone. 
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These problems highlight the need for careful analysis of the effects of noise and polariza- 
tion detection thresholds in studies of the fractional polarization of faint sources. Previous 
studies hav e focu se d on the fractiona l pola rization of bright radio sources in the NVSS. 
Mesa et all ^QM), IXucci et al.l J2004h . and ISeck fc Gaenslerl kooi ) considered polarized 
sources in the NVSS with Stokes / flux density stronger than 80 mJy, 100 mJy, and 80 mJy 
respectively. The high flux density thresholds in these studies were adopted to achieve a 
good level of completeness in 11 down to the limit set by residual instrumental polarization 
(n ~ 1%). These studies found a 11 distribution that decreases monotonically with increas- 
ing n, and with a median 11 ^ 1.8%. For these bright sources, noise effects are small, and the 
n distribution should be close to the intrinsic IIo distribution down to limits set by residual 
instrumental polarization. 

In this paper we investigate the shape of the IIo distribution for much fainter Stokes 
/ flux densities than those considered in previous work. The present data have angular 
resolution similar to the NVSS, so differences in the measured degree of polarization because 
of a difference in resolution are not expected. Our results can be compared directly with 
results based on the NVSS. 

Noise effects such as those illustrated in Figure H] were taken into account by a Monte 
Carlo analysis. A set of simulated catalogs was generated to accurately represent the effects 
of noise and the polarized flux density detection threshold in the data. Stok es / flux densities 



of sim ulated sources were drawn from the fit to observed source counts by IWindhorst et al. 



(Il990l ). We assume in these simulations that the Stokes / source counts of polarized sources 
have the same shape as those for all radio sources. This is a reasonabl e assumption be- 
cause ^ 80% of radio sources in the NVS S display significant polarization (IMesa et al.ll2002 



Tucci et al.ll2004l : iBeck fc Gaenslerll2004j ). The intrinsic Stokes I flux density is multiplied by 
the degree of polarization, IIq, drawn from an assumed IIq distribution to obtain the intrinsic 
polarized intensity Po- 

The intrinsic Stokes Jo and po of a simulated source are transformed to observed flux 
densities I and p by adding noise with statistical properties identical to the properties of 
the noise in the data. First, Po is converted into intrinsic Stokes Qo and Uo, assuming a 
random polarization angle. The error in the flux density is assumed to consist of a part that 
is proportional to the noise in the image at the location of the source, and a part that is 
proportional to the flux density of the source, added in quadrature. The error in the flux 
density S is evaluated as 

(3) 



where S represents the intrinsic flux density Jq, Qo, or Uo, cr is the rms noise in the image 
and Ci and C2 are constants. The value of Ci was determined from the rms variation 
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of th e flux density of bright sources in the ten fields after field registration (ITaylor et al. 
20031 ). Flux densities of sources brighter than 100 mJy varied by 2.5% (rms) over the ten 
fields. From this w e adopt C^ = 0.025. T he value of the constant C2 = 1.3 was taken 



to be that found by iRengelink et al.l (119971 ) from Monte Carlo simulations for flux density 
errors in the Westerbork Northern Sky Survey (WENSS). The value of a is different for each 
source, to represent variation of the noise with location in the mosaic. The distribution of 
polarized intensities of the resulting simulated sources, have the same statistical effects as 
the observed polarized intensities, including the effects of noise bias, the detection threshold 
and the variation of the noise with position in the mosaic. 

To derive the IIo distribution directly from the data, we compare the distribution of 
the data in a log(/) — log(p) diagram with the probability distribution for simulated source 
catalogues. Figure [5] shows the observed d ata points, and the model probability distribution 
assuming the Ho distribution function of iBeck fc Gaenslerl (120041 ). As expected, sources 
brighter than ~ 80 mJy are represented well by this model. However, fainter sources in our 
sample are more highly polarized than predicted by this distribution. This is clearly visible 
in FigureOfor sources with 10 < J < 30mJy, where an offset exists between the distribution 
of observed sources and t he predicted ri dge of maximum source density. A 2-dimensional 
Kolmogorov-Smirnov test (|Peacocklll983l ) rejected the hypothesis that the data were drawn 
from the simulated distribution at the 99.9% confidence level. 

The best fit IIo distribution was derived by fitting source probability density distribu- 
tions to the data in the log(/) — log(p) plane for trial IIq distributions. The trial IIq distri- 
butions were represented by a low -order Gauss-Hermite se r ies, a lso called a Gram-Charlier 
series, following the description of I Van der Marel fc Franxl (Il993l ). 



/(Ho) = exp 



N 



(4) 



We consider only modest deviations from a Gaussian, since previous work on bright NVSS 
sources suggests that the shape of the distribution is nearly Gaussian. Assuming that the 
Ho distribution peaks at zero and declines monotonically with increasing Ho, we use only 
the fourth term {i = 4) in Equation HJ which results in symmetric deviations. A coefficient 
/14 > means that the wings of the di stribution are stronger than th e wings of a Gaussian 
distribution, as shown graphically by IVan der Marel fc Franxl (119931 ). Higher order terms 
were not considered because of the limited size of our data sample at this time. The Ho 
distribution defined by Equation H] thus has two free parameters, the Gaussian standard 
deviation, au^, and the amplitude of the lowest order symmetric deviation from a Gaussian, 
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h^. The fits maximize the hkelihood L of the data as a function of these parameters, 



^ = n ^(^i'Pil^no,^4), (5) 

2=1 

with the probability of an observed for a given au^ and 

P{Ii,PiWn,,h4) = jj Y] e~^^ . (6) 

model , 

The product over i is over all sources in the data, A^data- Whereas the sum over j is over all 
simulated sources in the catalog, Mmodeb for a particular and h^. Typically Mmodci ~ 10^. 
The values aj^i and cTp^i are the observed errors in Stokes I and polarized flux density for the 
i*'^ source. 

The maximum likelihood Ho distribution was found through a grid search over the pa- 
rameter space. The best fitting model has parameters = 7%, and /i4 = 0.05. The 
probability density function of this model is shown along with the data in Figure El The un- 
certainty in the best-fitting parameters was evaluated empirically. Three hundred randomly 
selected samples, each containing on average the same number of sources as the observed 
sample, were drawn from the best fitting simulated catalog. Each of these samples was fit- 
ted with the maximum likelihood fit to evaluate the spread of the best-fitting parameters. 
Two thirds of these fits yielded a aji^ within 1% of the maximum likelihood value 7%. We 
conclude that au^ = (7.0 ± 1.0)%. The fitted value of is not independent of auo, because 
a larger /14 can partially compensate for a smaller ano- From the same 300 experimental fits, 
two thirds yielded a value < 0.1. Although the data are consistent with = 0, those fits 
with /14 contrained to be zero yield an average au^ = (8 ± 1)%. The data therefore suggest 
that the Ho distribution may be somewhat broader than a simple Gaussian with aji^ = 7.0%. 

The best-fitting model was also subjected to a 2-dimensional Kolmogorov-Smirnov 
test. The hypothesis that the data were drawn from the best-fit distribution in Figure M 
was rejected at t he 98 % confidence level. This is much better than the result for the 



Beck fc Gaenslerl (120041 ) distribution, but it is still suggestive that all the data are not well 
represented by the maximum likelihood model distribution. This is entirely the result of the 
fact that the best fitting Gauss-Hermite distribution does not fit the bright sources in the 
sample very well. We were able to produce a better fit by creating a set of hybrid simulated 



catalogs that use the iBeck fc Gaensleil (120041 ) distribution for brighter sources and our best 



fitting distribution for faint sources. The transition between these regimes was made smooth, 
with equal contributions from the two distributions at a flux density of 30 mJy. These hybrid 
catalogs fitted the data significantly better, with the best fitting model, using crn„ = 7.0% 
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and /14 = 0.05 at low flux density, passing the 2-dimensional Kolmogorov-Smirnov test. The 
median fractional polarization of the best fit IIo distribution for polarized sources with Stokes 
/ below 30mJy is (4.8 ± 0.7)%. 



3.3. Source Counts 

To derive source counts of the polarized sources, the effect of the varying noise level 
over the mosaic on source detection probability (completeness correction) was measured 
as a function of flux density by repeating the source detection on 1000 simulated images 
having the same noise characteristics and source density as the data. Separate simulations 
were performed for total intensity and for polarized intensity images. Each simulated image 
contained random Gaussian noise smoothed to the resolution of the DRAO images, and 
the same uniform rms amplitude as the data. For each simulated polarization image, two 
independent noise images were created to represent the Q and U images. Sources were 
pla ced at random positions , with flux densities drawn from the source-count curve derived 



by [Windhorst et al.l (Il990l ) between 0.1 and 500 mJy. Sources below the detection limit 
were included in the simulations to account for crowding in the field, and the possibility that 
faint sources are observed above the detection limit because of noise. The resulting synthetic 
Stokes / and polarized intensity images were then searched for sources in an identical manner 
as for the observed images. 

Figure [7] shows the derived differential source counts for total intensity (log - log 
/) and polarized flux density (log - log p) at 1420 MHz, normalized to the Euclidean 
expectation in the conventional way. The polarized source counts are also listed in Table [2J 
Counts were derived in bins starting at 500 /iJy to avoid sources with completeness correction 
greater than 10. This resulted in the removal of the fifteen faintest sources from the counts. 
The Stokes / source counts show good agreement with the Windhorst curve, although we 
find somewhat higher numbers around 10 mJy. The polarized source counts are nearly fiat 
in the flux density range observed, consistent with the increased fractional polarization of 
the faint radio sources. 

Figure [7] shows predicted polarized source count curves deriv ed by convolving Wind- 



horst's Stokes I counts with the strong-source Ho distribution from lBeck fc Gaenslerl (120041 ) 
(dashed curve), as well as our hybrid Ho distribution fit to the log /— logp di stribution in ^[X^ 
(solid curve). The observed polarized source counts show an excess over the Beck &: Gaenslerl 



( 120041 ) extrapolation for p < 3 mJy, and are consistent with the prediction based on the de- 
rived higher Ho for these sources. The data point at the lowest flux density lies marginally 
above our predicted curve. This may suggest a continuing trend toward even higher frac- 



- 12 - 



tional polarization in the sub-mJy population. The complete DRAO 30-field survey of the 
region, and other deep polarization surveys, will test this. 



3.4. Identification with Spitzer objects 

The positions of the polarized radio sources were examined in the Spitzer SWIRE images 
of ELAIS Nl. Accurate position s for the sources were obtained from the VLA FIRST 



survey images (IWhite et al.l 119971 ). which provides 1.4 GHz continuum Stokes / images at 
5" resolution with a 1-a sensitivity of 150 yuJy beam~^. The sensitivity of the FIRST images 
is sufficient to detect a Stokes I counterpart for every polarized source in the ELAIS Nl 
deep field. Approximately 35% of the polarized radio sources showed resolved structure on 
scales of 2" to 30" in the FIRST images. A polarized source may be associated with a radio 
lobe instead of the core of a radio galaxy. Visual inspection avoided misidentification in 
such cases. As a comparison to the polarized sources we also searched the SWIRE images 
for identification of sources with Stokes I fiux density larger than 1 mJy but no detectable 
polarization. 

In total, 54 polarized sources were unambiguously identified with Spitzer objects. An- 
other 17 polarized sources had uncertain identifications, i.e. there was more than one possible 
counterpart within the errors of the radio position. Two sources were not covered by SWIRE. 
The remaining 12% of the sources (10 objects) have no counterpart in the SWIRE images. 
Similar statistics resulted from the search for counterparts of the Stokes / sources with 
no detectable polarized emission. A wide range in fiux density and angular size was found 
among the identified Spitzer galaxies. Although some faint galaxies appear unresolved in the 
Spitzer images, the identified polarized sources seem to be mainly associated with elliptical 
galaxies. Three of the polarized sources have counterparts in the Spitzer images that are too 
faint to appear in the Spitzer ELAIS Nl source catalog. These sources were not included in 
the subsequent analysis. 

Figure [S] shows a near-infrared color-color diagram of Spitzer galaxies identified with 
radio sources that had catalogued fiux densities in all four IRAC bands at 3.6, 4.5, 5.8, and 
8.0 /zm. This includes 41 of the polarized sources. Since many of the sources are identified 
with extended galaxies, we used the isophot al fiuxes in each ban d. Our analysis of this color- 



color diagram is based on the modeling of ISajina et al.l (120051 ). which divided the diagram 



into the four regions separated by dashed lines in Figure O Boundaries between these 
regions were defined so as to separate galaxies depending on the strength of near-infrared 
PAH bands and the slope of the near-infrared continuum. Region 1 is mainly populated 
by sources with a continuum that rises with wavelength. The near-infrared spectrum of 
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these sources is usually dominated by non-equilibrium emission of stochastically heated very 
small dust grains, interpreted as PAH destruction by the hard ultraviolet spectrum of an 
AGN. Region 2 is occupied mainly by dusty star-forming galaxies with strong PAH bands 
at red-shift z < 0.5. This is because the strongest PAH features at low redshift contribute 
to the flux in IRAC bands at 3.6 and 8.0 fim. Region 3 is occupied by galaxies with fainter 
PAH emission, or by dusty starforming galaxies at redshift z = 0.5 — 1.5. A conspicuous 
concentration of galaxies in region 3 is the blue clump near log(S'5.8)/log(S'3.6) = —0.6, 
log(5'8.o)/ log(S'4.5) = —0.8. Galaxies in the blue clump have a near-infrared continuum that 
declines with wavelength, since the spectrum is dominated by the starlight of an old stellar 
population. These are elliptical galaxies, at a wide range of redshift. Region 4 is populated 
by PAH-dominated sources at redshift z = 1.5 — 2. 

The number of sources by region in Figure [S] is listed in Table [31 The host galaxies 
of the polarized sources occupy mainly regions 1 and 3. The source in region 4 cannot be 
considered a convincing high-redshift galaxy in view of uncertainties in the photometry and 
its proximity to the boundary with regions 1 and 3. The majority (71%) of the host galaxies 
of polarized sources in the Spitzer color-color diagram are either in region 1 or in the blue 
clump. In both cases, the radio emission is interpreted as emission from an AGN. 

Eleven polarized sources (27%) are found in the region of PAH-dominated galaxies, 
with log(S'8.o)/log(S'4.5) > —0.5, a few dex above the blue clump. Their location in the 
color-color diagram suggests either PAH- dominated galaxies at redshift 0.5 — 1.5, or galaxies 



with faint PAH bands at lower redshift (ISajina et al.ll2005l ). Inspection of the FIRST and 



SWIRE images of these sources shows that all appear to be elliptical galaxies with a smooth 
morphology, and some have resolved symmetric radio sources, suggestive of radio lobes. This 
suggests that the polarized radio emission in these galaxies is associated with AGN activity, 
while the PAH emission originates from dust at a substantial distance from the AGN, where 
it is shielded from the hard ultraviolet radi ation that would destr oy the PAHs. Dust in 



elliptical galaxies is a common phenomenon (iGoudfrooij et al.lll994l ). It can be produced in 



the envelopes of cool red giant stars or it can be acquired through a merger with a gas-rich 
galaxy. 

The Stokes / sources with no detectable polarized emission generally occupy the same 
areas of the color-color diagram as the polarized sources, with two exceptions. First, a 
significant fraction (15%) of the Stokes / sources is located in region 2, the area where star 
forming, PAH-dominated galaxies are expected. These objects are generally the faintest 
Stokes I radio sources, and none would be detectable in polarization at our sensitivity level. 
The Stokes / sources in region 2 are likely members of a population of star-forming galaxies 
that is believed to make up a large fraction of the radio source population below ~1 mJy. 
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The second exception is that Table [3] indicates an excess of polarized sources (27%) relative 
to sources with no detected polarization (14%) in region 3b, associated with galaxies having 
PAH emission. This difference is significant if Poisson errors are assumed. Confirmation will 
require a more complete sample of polarized radio sources identified with Spitzer galaxies. 



3.5. Nature of the mJy polarized source population 

The polarized sources found in the ELAIS Nl field have a median Stokes / flux density 
of 12 mJy. Models of radio source populations fitted to the total radio source counts suggest 
that most radio sources with 1420 MHz flux density > 1 mJy are steep-spectrum radio galax- 
ies, whose power is ultimately derived from accretion onto a compact object. However, it is 
not clear a priori that a faint polarization-selected sample of radio sources is representative 
of the entire population. In principle, a highly polarized population of faint radio sources 
may constitute a significant fraction of a sample of faint polarized sources. 

In Section [3721 we presented evidence that faint extragalactic polarized radio sources are 
on average more highly polarized than bright sources, with a median fractional polarization 
approximatel y three time s high er. From an analysis of sources with Stokes I > 100 mJy 



in the NVSS, iMesa et al.l (120021 ) also noted that the median fractional polarization of radio 
sources increases with decreasing flux density, from 1.05% for St okes / > 800 mJy to 1.84% 



between 100 — 200 mJy. From a similar analysis of the NVSS, iTucci et al.l (120041 ) showed 
that the anti-correlation between flux density and percentage polarization occurs only for 
steep-spectrum sources. They found that the median percentage polarization for steep- 
spectrum sources increased from 1.14% for flux densities greater than 800 mJy to 1.77% 
between 100 — 200 mJy. Our result extends this to much lower flux densities and indicates 
a much stronger effect for faint sources, resulting in a median polarization of 4.8% at Stokes 
/ = 10 — 30 mJy. Analysis of Tabled] shows that these faint polarized emitters are dominated 
by steep-spectrum sources; all but one of the polarized sources for which a 325 -1420 MHz 
spectral index exists (75%) has a < —0.4. 

Polarized sources that can be identified with galaxies in the Spitzer ELAIS Nl deep field, 
are associated with elliptical galaxies. Most of the host galaxies have near-infrared colors 
typical for dust emission from the vicinity of an AGN, or an old stellar population with 
no significant dust emission. The majority of the polarized sources is associated with AGN 
activity for this reason. The remaining 11 polarized sources associated with galaxies having 
PAH emission are also likely to contain an AGN. Some are clearly resolved double-lobed 
objects in the FIRST images, and all appear to be elliptical galaxies with a smooth brightness 
distribution. We conclude that there is no evidence for galaxies in our sample of polarized 
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sources in which the radio emission is powered by star formation. The higher fractional 
polarization of faint radio sources may be related to a population of radio-quiet AGN in 
which fainter radio emission correlates with conditions that favor increased polarization, 
for example more ordered magnetic fields or less internal Faraday depolarization. High- 
resolution polarimetry of these objects will provide more insight into their nature. 

4. Conclusions 

We present sensitive observations of a complete sample of compact polarized radio 
sources, as part of a deep integration of the ELAIS NX region made with the Synthesis 
Telescope at the Dominion Radio Astrophysical Observatory. A total of 83 polarized sources 
was detected in the ten-field mosaic. 

The distribution of fractional polarization of faint polarized sources was investigated 
with a Monte-Carlo analysis that generates synthetic source lists with the same noise statis- 
tics and observational selection criteria as the data. Maximum-likelihood fits of the synthetic 
source lists to the data in the log(/) - log(p) plane yielded a best fitting Gauss-Hermite func- 
tion with = (7.0 ± 1.0)%, /14 = 0.05 ± 0.05 for the distribution of intrinsic fractional 
polarization. The data demonstrate a trend of increasing fractional polarization with de- 
creasing flux density. 

Polarized source counts from the ELAIS Nl deep field are presented down to 0.5 mJy. 
We find that the Euclidean-normalized polarized counts remain flat below 1 mJy. The dis- 
tribution of fractional polarization derived from our Monte Carlo analysis is convolved with 
the total-intensity source counts to produce a prediction of the polarized source counts. The 
predicted Euclidean-normalized polarized counts are nearly flat to ~ 2 mJy, in good agree- 
ment with the data. However, the data at the faintest polarized flux densities suggests a 
continuing trend of increased polarization fraction with decreasing flux density. 

The near-infrared color-color diagram for host galaxies identifled with the polarized 
sources in the ELAIS Nl held shows that most of the host galaxies are eUipticals, or galaxies 
for which the near-infrared spectrum is dominated by stochastically heated very small grains, 
presumably from the vicinity of an AGN. Some of the host galaxies appear to have PAH 
bands in their near-infrared spectrum, but the morphological resemblance with ellipticals, 
and the fact that some of these polarized sources are resolved radio galaxies in the FIRST 
survey, indicates that these objects also harbour AGN. We suggest that the higher degree of 
polarization indicates a difference between AGN observed at a flux density of hundreds of 
mJy, and fainter AGN. 
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Right Ascension (J2000) 

Fig. 1. — Survey area of the first 10 fields of the DRAO mosaic in relation to the sky 
coverage of the original ELAIS Nl field (dashed lines) and the SWIRE survey (solid black 
lines). The ten field centers are indicated by the + symbols. Gray contours indicate the 
sensitivity at 1.5, 2, and 3 times the theoretical noise level in the center of the mosaic, and 
the edge of the field of view (thick fine) . 

Fig. 2. — [This Figure is provided as a separate image f2.gif] Continuum images of 
the ELAIS Nl field at 21 cm wavelength. Top row: Stokes I (gray scales linear from —0.1 to 
+5 mJy beam~^). Middle row: Stokes Q (gray scales linear from —1 to +1 mJy beam~^). 
Bottom row: Stokes U (gray scales linear from —1 to +1 mJy beam"^). Panels on the right 
show an enlargement of the area indicated by the white frame in the Stokes Q image. 
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Fig. 3. — [This Figure is provided as a separate image f3.gif] The distribution of 
amplitudes in the Q and U images after dividing by the mosaic weights to produce an image 
with uniform noise over the map equal to the noise value at the map center. Gaussian fits 
to the distributions (solid curves) were used to measure the map center rms at 78 //Jy. 
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Fig. 4. — Illustration of the effect of noise and the polarized flux density detection threshold 
on n as a function of flux density. This figure shows the variation of 11 with flux density 
for a simulated sample of sources, all 5% polarized, with random Gaussian noise added in /, 
Q, and U. The error bars are derived from standard error propagation, assuming the noise 
i s known. The curve show s the effect of polarization bias on 11, defined as a/ {Po + cr"^)/! 
(ISimmons fc StewartI 119851 ). 
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Fig. 5. — Distribution in log(/) - log(p) of observed sources ( open circles) , compared w ith the 
simulated distribution that assumes the Ho distribution of iBeck fc Gaenslerl (l2004l ). Gray 
scales and contours show the 2-dimensional probability density function of sources in the 
simulated catalog. The inner contours enclose 25% and 50% of the simulated sources, while 
the outer contour encloses 90%. Two lines mark the loci of sources with 11 = 1% (right) and 
n = 10% (left). 
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Fig. 6. — Same as Figure O now for the best fitting Gauss- Hermite Ho distribution, with 
an„ = 7.0% and h = 0.05. 
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Fig. 7. — Euclidean-normalized source counts for total flux density (triangles) and polar- 
ized flux density (circle s). The upper solid curve is the fit to observed source counts from 



Windhorst et al.l (Il990l ). The lower solid curve shows polarized source counts predicted by 



convolving the Stokes / source counts with the IT distribution derived from ELAIS Nl data 
for faint sources and the lBeck fc Gaenslerl (12004J ) distribution for bright sources, as explained 
in the text. The dotted curve shows polarized source counts de r ived by convolving only the 



Beck &: Gaenslerl ( 120041 ) distribution with the IWindhorst et al.l (Il990l ) source counts curve. 
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Fig. 8. — Spitzer near-infrared color-color diagram of ELAIS Nl radio source host galaxies. 
Stars: polarized sources; circles: sources with no detectable polari zed emission . The d ivision 
of this diagram into four regions and their interpretation follows ISajina et al.l (120051 ) and is 
explained in the text. 
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Table 1. ELAIS Nl Polarized Sources 



Source RA DEC I po PA Uq 0:325-1420 

Number (J2000) (mjy) (mjy beam-l) (°) (%) 
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± 


0. 


,11 


30 


2.0 


± 


0.2 


-0.79 ± 


0.03 


30 


16 


10 


27.12 




0.14 


54 


12 


54.0 


± 


1.6 


8.51 


± 


0.28 


1.24 




0. 


,11 


30 


14.4 


± 


1.3 


-0.43 ± 


0.17 


31 


16 


10 


57.84 


± 


0.08 


55 


35 


24.0 


± 


0.7 


19.46 


± 


0.53 


1.61 


± 


0. 


,11 


5 


9.7 


± 


0.7 


-0.63 ± 


0.06 


32 


16 


11 


00.48 


± 


0.06 


54 


42 


03.6 


± 


0.6 


29.41 


± 


0.76 


2.03 


± 


0. 


,11 


19 


7.8 


± 


0.5 


-0.74 ± 


0.04 


33 


16 


11 


20.40 


± 


0.08 


55 


28 


44.4 


± 


0.9 


18.22 


± 


0.50 


1.39 


± 


0, 


11 


3 


8.8 


± 


0.7 


-0.88 ± 


0.05 


34 


16 


11 


21.12 


± 


0.30 


54 


31 


55.2 


± 


3.3 


3.62 


± 


0.34 


0.43 


± 


0, 


10 


-5 


9.9 




2.4 


> -0.81 




35 


16 


11 


29.04 


± 


0.29 


55 


51 


36.0 


± 


2.6 


13.22 


± 


0.41 


0.44 


± 


0, 


10 


25 


4.3 


± 


1.0 


-0.87 ± 


0.06 


36 


16 


11 


37.92 


± 


0.29 


53 


59 


34.8 


± 


4.2 


13.82 


± 


0.44 


0.82 


± 


0. 


,10 


-18 


5.9 


± 


0.8 


-1.03 ± 


0.05 


37 


16 


11 


38.16 


± 


0.26 


55 


59 


52.8 


± 


2.7 


4.50 


± 


0.23 


0.48 


± 


0. 


,10 


14 


11.3 


± 


2.4 


-1.19 ± 


0.11 


38 


16 


11 


50.88 


± 


0.22 


55 


00 


54.0 


± 


1.0 


9.14 


± 


0.41 


0.82 


± 


0. 


,10 


14 


16.3 


± 


2.1 


-0.63 ± 


0.12 


39 


16 


12 


12.24 


± 


0.02 


55 


22 


48.0 


± 


0.2 


312.36 


± 


7.86 


20.37 


± 


0. 


,52 


29 


8.0 


± 


0.3 


-1.11 ± 


0.03 


40 


16 


12 


24.00 




0.13 


55 


26 


02.4 


± 


1.5 


9.57 


± 


0.34 


0.74 


± 


0, 


10 


32 


8.7 


± 


1.2 


-0.56 ± 


0.12 


41 


16 


12 


28.56 


± 


0.34 


55 


06 


46.8 


± 


2.0 


10.06 


± 


0.30 


0.37 


± 


0. 


10 


59 


4.1 


± 


1.1 


-1.13 ± 


0.06 


42 


16 


12 


31.68 


± 


0.31 


54 


18 


10.8 


± 


2.5 


9.81 


± 


0.36 


0.43 


± 


0. 


,10 


21 


4.5 


± 


1.1 


-0.73 ± 


0.10 


43 


16 


12 


35.28 


± 


0.04 


56 


28 


19.2 


± 


0.4 


176.88 


± 


4.51 


11.47 


± 


0. 


,31 


39 


7.8 


± 


0.3 


-0.96 ± 


0.03 


44 


16 


12 


47.52 


± 


0.33 


55 


02 


31.2 


± 


2.9 


7.68 


± 


0.25 


0.35 


± 


0. 


,10 


51 


4.9 


± 


1.4 


-0.74 ± 


0.12 
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Table 1 — Continued 



Source RA DEC / Po PA IIo 0:325-1420 

Number (J2000) (mjy) (mjy beam-l) (°) (%) 



45 


16 


12 


51.36 


± 


0.20 


56 


03 50.4 


± 


2.4 


2.59 


± 


0.19 


0.61 


± 


0. 


10 


25 


20.9 


± 


3. 


,7 


> -1.04 




46 


16 


13 


02.40 


± 


0.15 


54 


32 27.6 


± 


1.3 


6.48 


± 


0.25 


0.73 


± 


0. 


10 


56 


11.5 


± 


1. 


,7 


-0.86 ± 


0.12 


47 


16 


13 


16.80 


± 


0.36 


56 


08 13.2 


± 


3.8 


3.85 


± 


0.36 


0.56 


± 


0. 


,10 


-28 


14.1 


± 


2, 


,6 


> -0.77 




48 


16 


13 


19.20 


± 


0.09 


54 


16 40.8 


± 


1.0 


5.78 


± 


0.43 


1.34 


± 


0. 


11 


5 


31.5 




2, 


,8 


> -0.50 




49 


16 


13 


25.92 


± 


0.24 


55 


39 39.6 


± 


2.4 


3.06 


± 


0.18 


0.35 


± 


0. 


10 


17 


12.2 




3, 


,6 


> -0.93 




50 


16 


13 


26.64 


± 


0.10 


55 


15 46.8 


± 


1.1 


14.60 


± 


0.14 


0.80 


± 


0. 


,10 


14 


6.5 


± 


0, 


,9 


-1.19 ± 


0.04 


51 


16 


13 


28.80 


± 


0.27 


56 


17 49.2 


± 


2.1 


17.29 


± 


0.57 


0.87 


± 


0. 


,10 


26 


6.8 


± 


0. 


,8 


-0.53 ± 


0.07 


52 


16 


13 


30.72 


± 


0.06 


54 


27 21.6 


± 


0.6 


80.29 


± 


2.03 


2.67 


± 


0. 


,12 


15 


3.7 


± 


0. 


,2 


-0.85 ± 


0.03 


53 


16 


13 


36.72 


± 


0.36 


54 


11 16.8 


± 


2.7 


1.70 


± 


0.36 


0.59 


± 


0. 


,10 


1 


22.9 


± 


4, 


,2 


> -1.33 




54 


16 


13 


41.76 


± 


0.19 


56 


11 49.2 


± 


1.4 


101.05 


± 


7.62 


2.80 


± 


0, 


12 


22 


10.9 


± 


0, 


,6 


0.09 ± 


0.04 


55 


16 


13 


48.48 


± 


0.27 


54 


14 13.2 


± 


2.5 


8.12 


± 


0.40 


0.64 


± 


0. 


10 


18 


7.3 


± 


1, 


,2 


-1.10 ± 


0.07 


56 


16 


13 


56.16 


± 


0.25 


54 


57 28.8 


± 


2.5 


0.45 


± 


0.09 


0.28 


± 


0. 


,10 


-1 


50.6 ± 19. 


.2 


-3.97 ± 


0.03 


57 


16 


13 


56.84 


± 


0.22 


55 


02 08.2 


± 


2.3 


0.81 


± 


0.13 


0.50 


± 


0.02 


-31 


49.8 


± 5. 


,7 


> -1.82 




58 


16 


14 


00.96 


± 


0.20 


53 


57 21.6 


± 


1.8 


14.75 


± 


0.49 


1.08 


± 


0. 


,10 


32 


7.5 ± 0, 


,8 


-0.41 ± 


0.10 


59 


16 


14 


16.80 


± 


0.29 


55 


42 57.6 


± 


3.5 


1.14 


± 


0.36 


0.46 


± 


0, 


10 


56 


64.0 ± 17, 


,2 


> -1.60 




60 


16 


14 


21.12 


± 


0.06 


55 


36 39.6 


± 


0.6 


35.41 




0.93 


3.00 


± 


0, 


13 


18 


9.9 ± 0, 


,5 


-0.74 ± 


0.04 


61 


16 


14 


32.64 


± 


0.20 


55 


38 31.2 


± 


2.5 


1.27 


± 


0.39 


0.47 


± 


0, 


10 


-79 


90.8 ± 27, 


,5 


> -1.52 




oz 


ID 


iO 


z ( .00 


_|_ 
it 


U.lo 


54 


27 10.8 


_|_ 
it 


L.Z 




_|_ 
it 




u.y / 


_|_ 
it 


0. 


,10 




20.1 


± 


2, 


Q 
.0 


S> — U.DO 




63 


16 


15 


30.96 


± 


0.16 


54 


52 30.0 


± 


2.0 


5.43 


± 


0.22 


0.68 


± 


0. 


,10 


20 


13.1 


± 


2. 


.0 


-1.16 ± 


0.09 


64 


16 


15 


36.72 


± 


0.19 


53 46 37.2 


± 


2.4 


54.20 


± 


1.46 


2.64 


± 


0. 


,12 


19 


5.6 


± 


0. 


.3 


-0.85 ± 


0.03 


65 


16 


15 


49.68 


± 


0.06 


55 


16 40.8 


± 


0.7 


27.80 


± 


0.72 


1.75 


± 


0. 


11 


35 


7.4 


± 


0, 


,5 


-0.68 ± 


0.04 


66 


16 


16 


23.04 


± 


0.14 


55 


27 00.0 


± 


1.3 


13.10 


± 


0.40 


0.90 




0, 


10 


27 


8.8 




1, 


,0 


-0.89 ± 


0.06 


67 


16 


16 


23.52 




0.17 


54 


57 43.2 


± 


1.6 


10.35 


± 


0.48 


0.60 




0, 


10 


10 


6.0 




1, 


,0 


-0.97 ± 


0.07 


68 


16 


16 


37.92 


± 


0.18 


55 


45 14.4 


± 


2.1 


74.66 


± 


1.90 


0.85 




0, 


10 


51 


1.3 


± 


0, 


,2 


-0.61 ± 


0.03 


69 


16 


16 


39.36 


± 


0.11 


53 


58 12.0 


± 


0.9 


351.23 


± 


9.13 


5.14 


± 


0. 


,19 


23 


1.7 


± 


0, 


,1 


-0.85 ± 


0.03 


70 


16 


16 


40.08 


± 


0.18 


56 


20 38.4 


± 


1.5 


18.64 


± 


0.70 


1.37 


± 


0. 


,11 


24 


9.3 


± 


0. 


,8 


-0.52 ± 


0.07 


71 


16 


17 


57.60 


± 


0.22 


54 


51 36.0 


± 


3.0 


14.65 


± 


0.46 


0.68 


± 


0. 


,10 


-17 


5.9 


± 


0. 


.9 


-1.11 ± 


0.04 


72 


16 


18 


06.72 


± 


0.36 


54 


42 46.8 


± 


2.4 


5.31 


± 


0.26 


0.65 


± 


0. 


10 


24 


12.3 


± 


2, 


,0 


> -0.55 




73 


16 


18 


32.64 


± 


0.05 


54 


31 44.4 


± 


0.5 


48.68 


± 


1.27 


4.56 


± 


0, 


15 


-1 


10.3 




0, 


,4 


-0.78 ± 


0.03 


74 


16 


18 


57.57 


± 


0.12 


54 


29 26.2 


± 


1.3 


132.06 


± 


5.88 


2.55 


± 


0. 


,14 


21 


2.3 


± 


0, 


,2 


-1.39 ± 


0.03 


75 


16 


18 


59.28 


± 


0.18 


54 


52 40.8 


± 


1.7 


40.08 


± 


1.07 


1.86 


± 


0. 


,11 


21 


6.0 


± 


0, 


,4 


-0.50 ± 


0.04 


76 


16 


19 


15.36 


± 


0.20 


55 


05 13.2 


± 


1.4 


16.65 


± 


0.50 


0.69 


± 


0. 


,10 


-2 


5.7 


± 


0. 


,9 


-0.30 ± 


0.11 


77 


16 


19 


19.20 


± 


0.08 


55 


35 56.4 


± 


0.8 


53.48 


± 


1.36 


1.98 


± 


0. 


11 


15 


4.1 


± 


0, 


,3 


-0.44 ± 


0.04 


78 


16 


19 


19.44 


± 


0.06 


54 


48 25.2 


± 


0.6 


57.59 




1.46 


3.32 


± 


0, 


13 


11 


6.1 


± 


0, 


,3 


-0.67 ± 


0.03 


79 


16 


19 


24.24 


± 


0.18 


55 


50 52.8 


± 


1.6 


2.18 


± 


0.21 


0.79 


± 


0, 


10 


19 


26.8 


± 


3, 


,7 


> -1.16 




80 


16 


21 


13.68 


± 


0.10 


55 


23 42.0 


± 


0.9 


58.11 


± 


1.49 


2.62 


± 


0. 


12 


50 


4.8 


± 


0, 


,3 


-0.94 ± 


0.03 


81 


16 


21 


18.72 


± 


0.28 


55 


38 27.6 


± 


2.3 


7.97 


± 


0.40 


0.74 


± 


0. 


,10 


7 


8.5 


± 


1, 


,2 


-0.98 ± 


0.08 


82 


16 


21 


45.36 


± 


0.10 


55 


49 37.2 


± 


1.0 


31.10 


± 


1.01 


2.05 


± 


0. 


,11 


60 


6.9 


± 


0. 


,4 


-0.62 ± 


0.04 


83 


16 


22 


08.64 


± 


0.14 


55 


24 28.8 


± 


1.6 


14.23 


± 


2.15 


4.55 


± 


0. 


,15 


6 


27.7 


± 


1. 


,2 


-1.57 ± 


0.03 
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Table 2. Polarized source counts 



p 


N 


p^-^dN/dp 


(mJy) 




(Jyi-5 sr-1) 


0.71 


29 


1.38 ± 0.26 


1.42 


20 


1.26 ± 0.28 


2.50 


11 


1.10 ± 0.33 


4.61 


5 


1.01 ± 0.45 


9.94 


2 


1.46 ± 1.03 


20.4 


1 


2.32 ± 2.32 



Table 3. Radio sources in the IRAC color-color diagram 



Region 


Polarization detected 


Other sources 




Number 


% 


Number 


% 


1 


8 


20± 7 


32 


19±3 


2 








25 


15d=3 


3a^ 


21 


51±11 


80 


47±5 


3bb 


11 


27± 8 


24 


14±3 


4 


1 


2± 2 


9 


5±2 



^Region 3 blue clump, selected by 
log(5'8.o)/log(54.5)^-0.5 

''Region 3 PAH, selected by \og{Ss,Q) / \og{Si,^) > 
-0.5 
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